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Abstract — The lifting scheme of the discrete wavelet transform 
(DWT) is now quite well established as an efficient technique 
for image denoising. The lifting scheme factorization of 
biorthogonal filter banks is carried out with a linear-adaptive, 
delay free and faster decomposition arithmetic. This adaptive 
factorization is aimed to achieve a well transparent, more 
generalized, complexity free fast decomposition process in 
addition to preserve the features that an ordinary wavelet 
decomposition process offers. This work is targeted to get 
considerable reduction in computational complexity and power 
required for decomposition. The hard striking demerits of 
DWT structure viz., shift sensitivity and poor directionality 
had already been proven to be washed out with an emergence 
of dual tree complex wavelet (DT-CWT) structure. The well 
versed features of DT-CWT and robust lifting scheme are 
suitably combined to achieve an image denoising with prolific 
rise in computational speed and directionality, also with a 
desirable drop in computation time, power and complexity of 
algorithm compared to all other techniques. 

Index Terms — Lifting scheme, Dual Tree Complex Wavelet 
transform (DT CWT), Denoising, Computational complexity, 
coding gain and PSNR. 

I. Introduction 

One of the fundamental challenges in the field of image 
processing and computer vision is image denoising, where 
the underlying goal is to produce an estimate of the original 
image by suppressing noise from a noise contaminated 
version of the image. Image noise may be caused by different 
intrinsic (i.e., sensor) and extrinsic (i.e., environment) 
conditions which are Often not possible to avoid. Images are 
most frequently bear upon by noise evincing the facts like 
image capturing sensor internal imperfections, scarce of 
proper illumination of the object to be captured, during the 
process of its acquisition mean by digitization and its 
transmission. In fact the performance of the image acquiring 
sensors are generally affected by a wide variety of factors 
some among which are the factors such as environmental 
conditions during its acquisition process and the quality of 
the sensing elements deployed in image capturing systems, 
i.e., image acquisition sensors like Charge Coupled Device 
(CCD) cameras. For instance in acquiring images with a camera 
capturing object illumination and lighting levels and sensor 
temperature are major factors affecting the amount of noise 



in the resulting image. Another major possibility of image 
corruption with noise will be its transmission to a point of 
interest from a point of its perception by a sensor, principally 
due to the interference in the channel used for its 
transmission. Image noise is usually with reference to wide 
stochastic variations as opposed to deterministic distortions 
such as shading or lack of focus. It is actually the degree of 
variation of pixel values caused by the statistical nature of 
radioactive decay and detection processes. Even if we acquire 
an image of a uniform (flat) source on an ideal gamma camera 
with perfect uniformity and efficiency the number of counts 
detected in all pixels of the image will not be the same. In 
addition to noise added inherently by a sensor, image 
processing techniques also corrupt the image with noise [1]. 
There are a wide variety of noise models which can degrade 
the image, among which the most common is the additive 
one i.e., the additive white Gaussian noise will show a seri- 
ous impact on visual perception of image by its intensive 
degradation. And the remaining possible noise models which 
have a probability of corrupting an image are salt and pepper, 
speckle, poison and thermal noises in the communication 
medium. This undesirable corruption of image by noise is 
inherent in any image acquisition device. Certainly it is pos- 
sible to distinguish between two regimes when the degrada- 
tion structure was observed by photo sensors: in the first 
scenario, the measured intensities are sufficiently high and 
the noise is assumed to be signal-independent degradation. 
In the second scenario, only a few photons are detected, 
leading to a strong signal -dependent degradation. This 
image degradation by noise will make the image to be null 
information conveying object. Removing such noise is of 
great benefit in many applications and this may explain vast 
interest in this problem and its solution. The ultimate goal of 
image denoising technique is to improve the degraded image 
in some sense by suppressing the random noise, which has 
corrupted the image, while preserving the most important 
visual features of the image, such as edges. Denoising is an 
essential step prior to any higher-level image-processing 
tasks such as segmentation, photo restoration, visual track- 
ing where obtaining the original image content is crucial. 
Many algorithms have been proposed for image denoising, 
and there has been a fair amount of research on wavelet 
based image denoising, because wavelet provides an 
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appropriate basis for image denoising. But this single tree 
wavelet based image denoising has poor directionality, loss 
of phase information and shift sensitivity as its limitations. 
Hence we proposed dual tree wavelet structure for an im- 
proved quality image denoising while getting rid of all above 
limitations and lifting scheme structure for faster, complex 
free and efficient decomposition process. 

II. Dual Tree and Lifting Frame Work 

The conceptual idea central to a wavelet or sub-band 
transform [2], is illustrated in figure (1). The real wavelet 
transform makes an use of analysis filters h (n) and 
hj(n), followed by sub sampling, while the reverse transform 
first up-samples and then uses two synthesis filters f (n) and 
fj(n). Emulating same process with { g z Cry) , and 

{j%G*J hPiO*}} as analysis and synthesis filter pairs[3]with 
suitable sub sampling and up- sampling processes will lead 
to an imaginary tree wavelet transform, where the analysis 
and synthesis filter pairs in one tree will differ with those in 
another tree with a half sample delay. 

(7 D (vi"} ^ iJ D (u'}e _ ' |9! "- = w/2 (1) 

Since the filters in both trees can be made to offset by half 
sample, the wavelets resulting from the filter pair satisfies the 
Hilbert transform condition, given as 

iW r (w}; w > 0\ (2) 
r**V(w); w < O ) 
For real and imaginary trees, the modulation Matrices [1] 
indicated by M L (zl and M- (z~), respectively , defined as 

"--(kg &=3) 

The condition for perfect reconstruction can be represented 
in z-domain as 
For real 

H (-z)F (z) + /f 1 (-z)F 1 (z) = 
//„O)F U) + // 1 (z)F 1 (z) = 2d' 1 
For imaginary 

G (-z)P (z) + G 1 (-z)P 1 (z) = 
G (z)P (z) + G 1 (z)P 1 (z) = 2d' 1 
All the filters are assumed to be causal for simplicity and is 
the cause for introducing the term r" ; .secondly, (z) and 

M z < C- j are the dual versions of ."►/- - (z"l and m- z (zj respec- 
tively. We can write, for perfect reconstruction as 

M^Cz~ 1 ') t Af 10 GO = 2/ (7) 
M^z-^-y M 20 iz~) = 21 ® 



(5) 



(6) 




Fig 1(a): Real tree wavelet transform 
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Fig 1(b): Imaginary wavelets transform. 

Where I being an identity matrix of order '2' .If all filters are 
FIR, thenM 10 (s), JH^„ {z) and M 1L (z), M 21 (z) belongs to 
GL(2;R[z,£~*-]).As a special case in a orthogonal wavelet 
transform, F D (z) = H D (z), F L (z) = H L (z}, P D (z) = G c (z") 
and ?. (z) =t i (z), the modulation matrices 

M : Cz} = j i* * ( z \ and j - - \zj = jj-. (z), which are intern 2 
times a unitary matrix. In poly phase notation, these filters 
can be written interims of their even and odd phases, 
according to the following relations 

H„(z) = tfo^O 2 ) + z- 1 // 0o (z 2 ) 
/fi(z) = // le (z 2 ) + z^H^iz 2 ^ 
■■ G 0e (z 2 ) + z- 1 G 0o (z 2 
G le (z 2 ) + z^G^Cz 2 ) 



G (z) 
G x (z) 

Fi(z) 
P O) 
Pi 00 



Fo^z^ + z-^Cz 2 ) 



(9) 
(10) 
(11) 

(12) 



F le (z 2 ) +z- 1 F 1o (z 2 ) 

P le (z 2 ) + z- 1 P l0 (z 2 ) 
Where ff jj (zj i" contains even coefficients and H: (z~) con- 
tains odd coefficients for i=0,l.Thus 

fc fc 

«i»<>) = 2 ft i2i< z ^ anrf «i D <M = ^ fti 2JC+1 ^~ k (13) 

fc fc 

<^o e O) = ^^^o 2 ^~ k and G 0o <» = ^ go 2 A-+i z ~ fc 

GleC^) = ^l^ 1 ^^ = J]*!^-^" (14) 

fc fc 

fc fc 

*i.00 = 2 f*z**" nnd 'Sr-OO = ^fi X K+% z " (15) 

fc fc 

fc fc 

-Pi„GO = ^p l2fc z- fc anri P la GO = ^ Pi2 fc+ i^^ (16) 

fc fc 

We then assemble the polyphase matrices 



p ( \ — ro"( z J 
loC J _ U (z) /^ov 

So that 



(17) 
(18) 
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_ oGO[l _ Z J (19) 
^ 2o Cz 2 ) t = |iW 2oC z)[J (20) 

We now define ttTTz) and F^Tz"; as dual version of 

P ll? (z) and F :D Cz] ,then the perfect reconstruction property 
yields that 

P io (z)P ao (^l) t = I (2D 
P 20 (z)P 20 (z-^y = I (22) 
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Where P(z) and PTzt contain only Laurent polynomials. The 
resultant real and imaginary wavelet transforms in polyphase 
notation can be represented schematically in figure (2) 

<BH~ 





H" 1 h 


















Fig 2.a:Polyphase representation of real tree wavelet transform. 
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Fig 2b. Polyphase representation of imaginary tree wavelet 
transform. 

In both cases first the image has been sub sampled 
into its even and odd parts, then apply the dual polyphase 
matrices. For inverse transform case, first apply the polyphase 
matrix and then join the even and odd parts. The major prob- 
lem of an FIR wavelet transform resides in finding the matri- 
ces r-r(jjand ?--(-) with determinants as one. Once we 

have such a matrices, >C7Tz'\ and p^7Tz"\ the four filters for 
the wavelet transform will be derived immediately from 
polyphase matrices and . From Cramer's rule 

// 0e (z) = F lo (z" 1 ) , H 0o (z) = -F le {z _1 ) 



= F 0o (z J ) and H lo (z) = he^ 1 ) 
G 0e( z ) = ^10 O" 1 ) and G 0o( z ) = -Plei 2 ' 1 ) 

G te (z) - -PqoI^ 1 ) and G io( z ) = Pa^ 1 ) 
This implies 

H^z) = z^-z- 1 ) and H (z) = -r ^(-z" 1 ) 
G t (z) = z-%{-z x )ani G (z) = -z^P^-z" 1 ) 



(23) 
(24) 

(25) 
(26) 



are more general in the sense that all the classical wavelets 
can be generated by the lifting scheme. For faster 
implementation of wavelet transform,the lifting scheme makes 
an optimal use of the similarities between the high pass and 
the low pass filters. The flexibility afforded by the lifting 
scheme allows the basis functions associated with wavelet 
coefficients near the window' boundaries to change their 
general shape .This new basis function can be used to 
minimize boundary effects. Being a new method for 
biorthogonal wavelet construction, the lifting scheme allows 
the faster implementation of wavelet transform. The number 
of flops can be reduced by a factor of two. The basic block 
diagram of lifting steps is illustrated in figure .3 




Fig. 3: Representation of lifting steps. 

Where P r is a predict [7] operator and U is an update 
operator [8]. The lifting scheme can also be used in situa- 
tions, such as wavelets on bounded domains, wavelets on 
curves and surfaces, weighted wavelets, and wavelets with 
irregular sampling. Conceptually the lifting scheme originates 
with a 'lazy wavelet', a function which essentially doesn't do 
a thing , but having a formal properties of a wavelet. The 
lifting scheme then generally builds a new wavelet, with im- 
proved properties by adding in, a several new basis func- 
tions to the lazy wavelet to make it faster as shown in fig- 
ure(4). This is the inspiration behind the name 'Lifting 
Scheme' . 

'il' 



The most trivial example of a polyphase matrix is P(z)=I. 
This will result in, F~ (zj = H e (z~) = land 

F L (z) = H. (z) = z~- ■ Thus the wavelet transform does 
nothing else but sub sampling even and odd samples. This 
transform is generally called as a polyphase transform, but in 
the context of lifting it is often referred to as the lazy wavelet 
transform. This lazy wavelet transform [4] incorporates a 
numerous complex computations, which will consume a large 
amount of power, and time for computation, which results in 
high coding and computational complexity. 

To avoid these hard striking shortcomings due to process 
complexity and delay, the natural wavelets are further 
configured to yield a simple, flexible and computationally 
efficient second generation wavelet structure.. The robust 
lifting scheme [4] [5] was introduced to build second 
generation wavelets referred commonly as lifted wavelets, 
which are not necessarily translates and dilates of one fixed 
function [6]. This construction is entirely spatial and therefore 
ideally suited for building second generation wavelets, that 
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Fig. 4: Fast wavelet transform using lifting scheme 

which first does a lazy wavelet transform and then computes 
Yj-Urr*. and finally lifts the wavelet coefficients of the input 
image 

The lifting scheme uses the relationship between the 
perfect reconstruction filter pairs {/ : (»J . f, (n)J,that have the 
same low pass or high pass filter. The filter pairs 
{&faXit£flQ } an d {p^(f0ipi (")} to be complementary if in 
case the corresponding polyphase matrices P 1B (z) and P 2 p (z) 
have their determinant as 1 . and if the filter pair, { } and { } are 
complementary, then consequently the filter pairs { } and 
{ }are also complementary. Then any other finite filters and, 
complementary to and respectively are of the form 
'')"""'(/) = FiO) + f D (z)5(z 2 ) (27) 
Pi new (z) = Pi GO + PqGO^ 2 ) (28) 
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Where S(z) is a Laurent polynomial. Similarly as per as dual 
lifting concerned, there exists any finite filters of the form, 

F a "'^Cx) = F QO + F t G0t(> 2 ) (29) 



r-i new 



GO = F GO + F a GOt(2 2 ) (30) 



Where t(z) is a Laurent polynomial. The resulted lifted 
wavelet tree structure for the case of classical sub band 
scheme ,when lifted the lowpass sub band with the help of 
the high pass sub band for both real and imaginary trees is 
shown in figure(5).And in a similar way the dual lifted wave- 
let tree structure for the case of classical sub band scheme, 
when lifted the high pass sub band with the help of the low 
pass sub band for both real and imaginary trees is shown in 
figure(6). 
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Fig .5. The Lifting scheme (a) real (b) imaginary trees: first a 
classical sub band filters scheme and then lifting the low pass sub 
band with the help of the high pass sub band. 




Fig. 6:The dual lifting scheme (a) real,(b) imaginary wavelet trees: 
First a classical sub band filter scheme and later lifting the high pass 
sub band with the help of the low pass sub band. 

III. Proposed Work 

In this paper, it is proposed to reduce the computational 
complexity, computational time, power and complexity of 
coding structure, with a substantial growth in average 
computational speed, and coding gain in a traditional dual 
tree complex wavelet transform, which employs two separate 
Quadrature distinct discrete wavelet structures in which one 
is taken as the real tree wavelet structure[9] [ 10], and the other 
is taken as the imaginary tree wavelet structure. Each discrete 
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wavelet transform structure irrespective of its nature will have 
a pair of low pass and high pass analysis filters and a similar 
pair of low pass and high pass synthesis filters. The analysis 
and synthesis filter pair in one tree will differ with those in 
other tree by an half sample delay. But the wavelets in the 
traditional system are lazy wavelets(i.e., computationally 
complex, slower and consumes a high quantity of 
power). Hence to modernize the wavelets and to cater the 
hard striking inabilities of traditional wavelets. The traditional 
wavelets are further processed and modified, to lift them, in a 
view to get a more simple, flexible, linear, complexity free and 
fast second generation wavelets[l 1]. These second 
generation wavelets are constructed, first by lifting the lazy 
wavelet which involves bi orthogonal factorization [12] [13] 
[14] [15] of lazy wavelet filters to eradicate delay artifacts and 
to remove the coefficients in the filter transfer functions 
responsible for delayed processing, or operational complexity, 
consuming higher processing power, to yield the lifted filters 
for second generation wavelets, called lifted wavelet filters 
for computationally faster lifted wavelet transform. 

The process of the lifting wavelet construction involves 
factorization of wavelet filters to get simple lifted wavelet 
filters. This activity involves running Euclidean algorithm 
starting from polyphase components of F- (z) and F z (z~) as 

{F-o iz~) > F Ds (z)} and { F f (z] , F Do (z"i} respectively and their 
gcd's are monomials. With the non-uniqueness of the divi- 
sion, we can always choose the quotients q.-(z),so that the 
gcd is a constant. Let this constant be then we thus have 

12.13 nri" oii'o-i (3D 

r = i 

15:81 nr"^ sra (32) 

in case jF Do (z}| > \F 0e {zj ,the first quotient q. (t Vis zero. 
We can always assume that 'n' is even. If we know the filters 
F P (z) and P c (z) ,then we can always find the complementary 

filters?; : (z) and F 1 P (z) by letting 



p 10 "00 



F00O) 



•□.GO PiAz)]rTr<ft( 

■coo p 10 °oojy L 1 
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K n 



(33) 



(34) 



Here the final diagonal matrix follows from the fact that 
det (p, Q z (z"\ jj=det(F-~ " (zl)=l and 'n' is even, then by rewrit- 
ing the above equations, we observe that 



\ qi [z) 11 = rl q,-(z)j ro li 



1 



(35) 



Using the first equation of above in case 'i' is odd and the 
second in case is even yields 





i 



(36) 



As an example to this factorization process to built 
computationally speed lifting wavelets, let us consider a lazy 
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wavelet of haar type. Now our task is to factor the filters 
associated with lazy haar wavelet transform to build , active 
lifting wavelet filters there by to build lifted haar wavelets. In 
the case of unnormalized wavelets we know that for a single 
tree 

F O) = 1 + 2- 1 , F % {z) = =± + \z- % , 

using the aforementioned Euclidean algorithm, we can thus 
write the polyphase matrices as 

In 

ll 



.OO = 



,c» = 



?] [o 



2 



f] 



2 



Thus on the analysis side we have 

P 10 (zT 1 = PiJ-) = \ 1 2 
10 1. 



«o--MI)-G TIL 1 , a 



(37) 



(38) 



(39) 



(40) 



Thus the forward and reverse wavelet transform 
implementation structure using the alternating lifting and dual 
lifting steps are better illustrated in figure.7 




Fig .7a .Forward wavelet transform using the lifting scheme. 
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Fig.7.b:The inverse wavelet transform using the lifting scheme 

d t = d t — Si Si ~ S, 
Similarly the inverse factorization yields, 

cZ^ — cij ~\~ . 



dl — ^27+1 
1 

2 



(41) 



- - di 

2 ' 




P 

•d' 

'P' 
•d' 

[ 



[ 0.27351197 0.27351197] 
[ 3.87782215-3.26868661] 
[ -0.28650326 -0.28650326] 
[-0.54859417 2.94176754] 



[01 
[2] 
[-2] 
[4] 



[ 0.09982322-0.34381326-0.34381326 0.09982322] 

[-2] 



[ 0.86857870] [ 1.15130615] 



[73] 



The transfer functions of the corresponding lifted wavelet 
analysis filters for the dual tree complex wavelet transform 
are given in figure (8) 



o.s 




10 
G1(n) 



Fig .8:Transfer functions of lifted wavelet dual tree analysis filters. 

Similarly the transfer functions of the corresponding lifted 
wavelet synthesis filters for the dual tree complex wavelet 
transform are given in figure .9 




X2l+1 =d l ° , x 2; = S ; ° (42) 

Thus the lifting wavelets are designed for both real and 
imaginary trees ,and incorporated to get real and imaginary 
discrete wavelet transform [1 1] trees, which are suitably 
combined to form the dual tree complex wavelet transform to 
erase the limitations of the single tree discrete wavelet 
transform, such as shift variance, poor directionality, loss of 
phase information and aliasing. 

Thus for the case of a boir 6.8 wavelet type, the corresponding 
lifting scheme LS = { ... 
'p' [-2.65899636 0.99715069] [0] 
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Fig. 9:Transfer functions of lifted wavelet dual tree synthesis 
filters. 

The corresponding lifted wavelet coefficients are given 
for the dual tree complex lifting wavelet transform both for 
analysis and synthesis separately as in table(l) and table (2). 

Thus a computationally fast dual tree complex wavelet 
transform is constructed based on lifted wavelets .Now this 
efficient and faster DTCWT tool is employed as a tool of 
interest for an image denoising. In this paper a noise degraded 
image is first applied to both real and imaginary discrete 
wavelet transform trees to obtain real and imaginary 
coefficients. These coefficients are suitably combined to get 
complex wavelet coefficients. 

The corresponding dual tree wavelet structure is shown 
in figure. 10 
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Table I. Analysis lifted wavelet coefficients for DTCWT with LS. 
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Table II. Synthesis lifted wavelet coefficients for DTCWT with LS 
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2D Dual-Tree Complex Wavelets 




Fig.l0:Dual tree complex lifted wavelet structure. 



The dual tree complex wavelet transform embedded with 
lifting scheme is initially aimed to achieve a noteworthy 
progress in the computational speed up in the image analy- 
sis. This computationally outstanding tool integrated with 
the globally accepted adapted wavelet threshoulding 
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techniques have been employed as the key techniques to 
develop a modified denoising algorithm for image denoising. 
Since the dual tree complex wavelet transform with lifting 
structure incredibly accelerate the decomposition process 
with the conserved shift invariance and directional selective 
properties as earlier , adaptive soft thresholding technique 
central to the adaptive wavelet threshoulding has been inte- 
grated to denoise the image. The algorithm is intelligibly pre- 
sented and implemented as 

1. Select the test image of size 256X256. 

2. Select the type of the wavelet to be lifted. 

3. Lift the selected wavelet with sufficient predict and update 
operations to generate the lifted wavelet structure of the 
selected wavelet . 

4. Design and generate the lifted wavelet filters associated 
with the lifted wavelet(Both for analysis and synthesis ). 

5. Apply the dual tree based complex wavelet transform on 
the test image using the lifted wavelet filter pair for 
analysis. 

6. Grouping the connected components using the 
neighborhood based statistical modeling of the complex 
wavelet coefficients. 

7. Apply the adaptive soft thresholding technique with a 
threshold value of 40. 

8. Applly the inverse dual tree complex wavelet transform 
on the thresholded complex wavelet coefficients using 
the lifted wavelet filters for synthesis. 

9. Compare the results obtained in terms of PSNR and RMSE. 

A. Computational complexity 

In this section ,we used to focus on the computational 
complexity of the wavelet transform .As a comparison base 
we used the standard algorithm(DTCWT), which 
corresponds to applying the polyphase matrix. This will take 
the advantage of the fact that the filters will be sub sampled 
and thus avoids computing samples that will be sub sampled 
immediately. For each filter either it may be low pass or high 
pass,the number of multiplications and additions required 
are \ +■ land ii? D (z}l respectively. Thus in a standard 

algorithm without applying lifting scheme, the number of 
multiplications and additions required for a single discrete 
wavelet transform tree are (!/f D (z}| -{-|fi L (z)l + 2} and 
(tlfi[}(z}l + l^iKl 13 respectively. Similarly for an imaginary 
tree, the ( G^f z) ! + 'G. (z) +■ 2)number of multiplications 
and( G- (z"\ — G- {:) ) number of additions are required. 
Thus for dual tree complex wavelet transform, the total number 
of multiplications and additions needed 
are(lff (z)| + \H t (z)\ + |^ (z)| + |C L (z)| + 4-) and 

(lfl D (z)l + 1^0)1 + IC D (z)l + iGiOOD respectively. For 
the case of symmetric filters let 
\H (z)\ = 2L t = \G fz)| , Ifl^z)! = 2L 2 = \G L (z}i , a n d 
also be assumed that L 2 > I 1 .Then for a standard algorithm 
,the number of multiplications and additions required for a 
real tree are 2l-+2l-+2 and 2L, +■ 2 L-. .Similarly for an 
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imaginary tree the total number of multiplications and 
additions required are 2+2+2 and 2. Then for the dual tree 
complex wavelet transform, the total number of multiplications 
and additions required are 4(4 and 4() respectively. Similarly 
if the lifting scheme is incorporated, then for a single tree, 2 
scaling operations are required, for lifting steps,4 lifting 
operations are required and for final lifting steps, the number 
of operations required are 2(+l).Then the total number of 
operations required are 2(), among which the total number of 
multiplications and additions are and respectively. Thus for 
the dual tree complex wavelet transform using the lifting 
scheme the total number of 2() multiplications and 2(additions 
are required. Correspondingly the computational operations 



required for various transform cases are investigated in 
table. 3. 

Graphically the computational complexity of the proposed 
algorithm with DTCWT-LS is compared with DTCWT as 
shown in fig. 1 1 . As computational complexity increases, the 
computational time increases and vice-versa. Thus an algo- 
rithm with higher computational complexity will take higher 
time for its computation and that with lower computational 
complexity will take less time for its computation on the same 
processor. Thus the computational time with several wave- 
lets are compared in terms of no of computations required as 
in table .4 



Table III. Comparison of computational complexity between standard algorithm and lifting algorithm. 



S NO 


Transform o r Tool 


No. of multiplications 


No of additions 


Total number of operations 


1 


Real DWT 


2(1, + i : ) + 2 


2(L-„ 


4(1, +■ £»> + 2 


2 


Imag DWT 


2(1, + L.-y + 2 


2(1, + I.) 


4(1, + + 2 


3 


Dual CWT 


+ I.) + 4 


4(1, - 


8(1, +■ L 2 J - 4 


A 


Real DWT with LS 


r,, j_ _ 2 


*» " ^ 


2(1, +• £-) - 2 


5 


Imag DWT with LS 


+ L z - - 




2(1, + L z ) + 2 


6 


Dual CWT with LS 


2(1, i 2) 


2(1, - 1;) 


4(1, - 4 



CODrXO PERFORMANCE IXTERMi Or COiP-.TA—OyjC COMPLEXITY 




TOOL ...ITKCT: TYPE 



Fig. 11 

Table IV Comparison of computational time between standard algorithm and lifting algorithm for various wavelet filters 



Wavelet 


1 1 I < V\ I 


DTCWT LS 


Growth (<?.'<.) 


1 1 : : : : : 


3 


3 


A> 




14 


9 


56°/o 




23 


14 


57% 


(9-7) 


23 


14 


64% 


( 4, 2)B spline 


17 


io 


70% 


Int erp o 1 at ing 


3(1, — 


3/2(Ll+L2) 


^= 100% 




4(1, +■ L z j +■ 2 


2(1, — li+2) 


5= 100% 



B. Computational Time: 

As the number of operations are increases then the 
computational time will also increases. Thus the 
computational time of a standard and lifting algorithms are 
compared as in table. 5, .Inversely the computational time of 
the dual tree complex wavelet transform with lifting scheme 
will be compared graphically with that of DTCWTas shown 
in figure. 12. 

C. Coding gain: 

As the number of operations in the lifting scheme based 
real imaginary and dual tree complex wavelet transform are 
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Table V 



S_No 


Tool 


Total computations 


Time 


1 


Real DWT 


4(1, - L.i - 2 


100% 


2 


Imag DWT 


4(1, + 1-j -f 2 


100% 


3 


Dual CWT 


8 Cl; - ! = > + * 


200% 


4 


Real DWT with LS 


2(L, +.£-•) h- 2 


50% 


5 


Imag DWT with LS 


2(1, +- 1.) + 2 


50% 


6 


Dual CWT with LS 


4(1, - L a : - 4 


100% 



exactly half of those in standard real, imaginary and dual tree 
complex wavelet transforms respectively, then the overall 
coding gain will increases by 100%. Graphically the coding 
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gain of the DTCWT-LS is compared with that of the DTCWT 
as shown in figure 1 3 

CODING PERFORMANCE IN TERMS OF EXECUTION TIME 



I 1.5- 



IA-.G TREfc W/0 



- EXECUTION TIME WITHOUT LS 
EXECUTION TIME WITH LS 



1 



1.2 1.4 IS 1.8 2 22 24 2.6 2E 
TOOL/METHOD TYPE 



Fig(l 2):Compaison of computational time of DTCWT-LS with 
that of the DTCWT 
CODING PERFORMANCE IN TERMS OF CODING EFFICIENCY OR GAIN 

35 



S 25 - 



it 15 

LU 

a 



DUAL WITH 



IMAG WITH 



- CODING EFFICIENCY W/0 LS 

- CODING EFFICIENCY WITH LS 



DUAL TREE W/C 



IMAG TREE W/O 




algorithms have been compared with several images of dif- 
ferent size and format and the results are given in table.7 .For 
instance with the DTCWT without LS algorithm the no of 
complex computations required to process a cameraman im- 
age of size 256X256 and of tif format, a real tree requires 
9568256 computations ,for an imaginary tree 9568256 com- 
putations are needed and for dual tree on a whole 19136512 
computations are essential. But with DTCWT with LS algo- 
rithm ,to process the same image it requires 4849664 real, 
4849664 imaginary and 9699328 dual tree computations are 
sufficient. Similarly to process a Barbara image of jpg format 
and of 256X256 size the DTCWT without LS algorithms takes 
a time of 0. 1486 seconds for real tree, 0.024 seconds for imagi- 
nary tree and on a whole 0.2106 seconds are essential for 
dual tree. Where as to process the same image the DTCWT 
with LS consumes a time of 0. 1239 seconds for real tree,0.0188 
seconds for imaginary tree and on a total 0. 1436 seconds for 
dual tree is enough. The testing results have strongly wit- 
nessed the computational dominance of the DTCWT with 
LS algorithm over the DTCWT without LS algorithm, and 
outstandingly judged the computational efficiency of the 
DTCWT with LS algorithm In order to compare the quality 
of the denoised image signals obtained with both the algo- 
rithms such as DTCWT without lifting scheme and DTCWT 
with lifting scheme, the percent root mean square formula for 
both original and reconstructed image signals is 



Figure 13 

IV. Results and Discussion 

The denoising results are presented and compared with 
the former algorithms employing the conventional DTCWT 
and adaptive wavelet thresholding. The simulation results 
obtained have ascertained that there is a remarkable progress 
in the PSNR value and notable drop in the RMSE value of the 
algorithm compared to the former algorithm. Practically a 
PSNR value of 40.04 and an RMSE value of 05.0267 are pos- 
sible for additive white Gaussian noise(AWGN),but with the 
former algorithm the maximum PSNR value of 22.8709 and 
minimum RMSE value of 14.9928 are only possible with 
AWGN. With an impulse noise ,the modified algorithm offers 
a maximum PSNR of 47.4560 and minimum RMSE of 
04.5656, where as the former algorithm can only offer a maxi- 
mum PSNR of 24.2356 and a minimum RMSE of 12.5612.. An 
important note here is that all the results are obtained with a 
Barbara.jpg image of 256X256 as a test image. 

A comparative analysis of the performance offered by 
the DTCWT without lifting scheme(LS) and the DTCWT with 
lifting scheme are done in terms of the no of computations 
required and the amount of time elapsed to perform them for 
a particular image format. The performance of the both the 
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PRD = - ■ X100 (43) 

Where x, is the denoised image signal obtained with the 

DTCWT without LS and x i is the resulting denoised image 
signal obtained by the DTCWT. Similarly the outputs ob- 
tained with the DTCWT with LS and DTCWT without LS are 
here subtracted from the original input and differences are 
computed as RMS error (RMSE). The generalized formula for 
RMSE is 



RMSE = 



(44) 



where x\ is denoised image signal and x, is the original 
image signal without any noise. 




(a) (b) (c) 

Fig.l4:Visual Results obtained with Gaussian noise and (a). Original 
image(b)Noisy image (c)Denoised image. 

The denoising results obtained with Gaussian noise and 
impulse noise are illustrated in figure. 14 & 15. 

Image quality performance metrics PSNR and MSE are 
computed and are tabulated as given in below table. 6. 
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(b) (c) 
Fig. 15 Visual Results obtained with an Impulse noise (a). Original 
image, (b)Noisy image (c)Denoised image. 

Table VI 



Noise Type 


Tech nique 


PSNR 


RMSE 


AWGN 


DTCWT 


32.0637 


6.3585 


DTCWT-LS 


47. 8222 


15.5608 


Impulse 


DTCWT 


22. 8946 


18.2731 


DTCWT-LS 


48.4823 


12.9214 



V. Conclusions 

In this paper we developed a algorithm which suitably 
combines the salient features of nearly shift invariant, 
directionally selective dyadic decomposition tree based dual 
tree complex wavelet transforms deploying quadrature dis- 
tinct ,perfect reconstruction filter banks and linear adaptive, 
and faster decomposition arithmetic based robust lifting 
scheme for an image denoising in a view to get a substantial 
speed up in the decomposition and reconstruction structure 
of the dual tree complex wavelet transform while preserving 
all performance salient features that results from the dual tree 
complex wavelet transform itself. This algorithm provides a 
substantial drop in the average decomposition time, average 
power required for decomposition as a major distinction fac- 
tor from an ordinary dual tree complex wavelet decomposi- 
tion algorithm. With this lifted Dual tree structure we achieved 
an accountable growth in the coding efficiency and coding 
gain. This algorithm first performs the factorization of the 
filters associated with the selected wavelet type based on 
the Euclidean algorithm and then constructs its associated 
lifting scheme from which the filters associated with the lifted 
wavelet will be derived. These filters are employed to achieve 
an in-place-implementation of the fast Dual tree Wavelet 
structure. This computationally efficient, dual tree complex 
wavelet structure has been employed as a tool of interest for 
an image denoising and achieved a considerable growth in 
image quality and performance metrics viz. , PSNR and RMSE. 
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Table VII Tested Results for Computational Complexity 





P-fa of Operations 


Computation time 


Image 


Dual TreeCWT 


For Lifting scheme DTCWT 


Dual Tree CWT 


DTCWT-LS 




Real 


Imaginary 


Dual 


Real 


Imaginary 


Dual 


Real 


Imaginary 


Dual 


Real 


Imaginary 


Dual 


cameraman 


9568256 


9568256 


19136512 


4849664 


4849664 


9699328 


0.547 


0.547 


0547 


0.1293 


0.0202 


0.1562 


peppers 


28704768 


28704768 


57409536 


14548992 


14548992 


29097984 


0.1712 


0.044 


0.3051 


0.1051 


0.0309 


0.1203 


Lena 


5840000 


584G000 


11680000 


2960000 


2960000 


5920000 


0.1455 


0.0194 


0.1951 


0.1212 


0.0162 


0.1629 


Brandyrose 


8380400 


8380100 


16760800 


4247600 


4247600 


8495200 


0.1515 


0.0217 


0.2073 


0.1241 


0.0175 


0.1402 


Barbara 


9568256 


9568256 


19136512 


4849664 


4849664 


9699328 


0.1486 


0.024 


0.2106 


0.1239 


0.0188 


0.1436 


Bear 


8906000 


8906000 


17812000 


4514000 


4514000 


9028000 


0.1536 


0.0229 


0212 


0.1233 


0.0182 


0.1313 


Kid 


8438800 


8438800 


16877600 


4277200 


4277200 


8554400 


0.1455 


0.0219 


0.2047 


0.1229 


0.0177 


0.1763 


Watch 


4380000 


438G000 


8760000 


2220000 


2220000 


4440000 


0.1454 


0.0189 


0.1928 


0.1205 


0.0154 


0.1564 


Parrot 


9568256 


9568256 


19136512 


4849664 


4849664 


9699328 


0.1479 


0.0236 


0.2088 


0.125 


0.0201 


0.1451 


Baboon 


5840000 


584G000 


11680000 


2960000 


2960000 


5920000 


0.1442 


0.0201 


0.1942 


0.103 


0.0162 


0.1192 


House 


9568256 


9568256 


19136512 


4849664 


4849664 


9699328 


0.1582 


0.0247 


0.2199 


0.1254 


0.0196 


0.1396 
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